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We introduce a new numerical technique, namely triangular tessellation, to calculate the free 
energy associated with the adsorption of a colloidal particle at a flat interface. The theory and 
numerical scheme presented here are sufficiently general to handle non-convex patchy colloids with 
arbitrary surface patterns characterized by a wetting angle, e.g., amphiphilicity. We ignore interfacial 
deformation due to capillary, electrostatic, or gravitational forces, but the method can be extended 
to take such effects into account. It is verified that the numerical method presented is accurate and 
sufficiently stable to be applied to more general situations than presented in this paper. The merits 
of the tessellation method prove to outweigh those of traditionally used semi-analytic approaches, 
especially when it comes to generality and applicability. 

PACS numbers: 82.70.Dd, 02.60.Gf, 68.03.Cd 



I. INTRODUCTION 

Small particles at liquid-liquid interfaces arc of scien- 
tific interest, but can also be exploited for industrial ap- 
plications. Particles adsorbed at an interface have a mul- 
titude of applications, ranging from the formation of two- 
dimensional structures [JJ [2] , which may be utilized in op- 
tical devices, to the stabilization of foams and pickering 
emulsions [3]. The range of sizes, shapes, and material 
properties with which colloids can be endowed, makes 
them the ideal constituents for self-assembled macro- 
scopic structures. In addition, colloid tunability allows 
tailoring to specific systems, which gives tremendous ad- 
vantages over atomic materials. A more fundamental im- 
petus to the study of colloid adsorption, is based on gain- 
ing a better understanding in phase transitions and criti- 
cal phenomena of two-dimensional fluids of nanoparticles 
at an interface. 

Many theoretical investigations of colloids at an inter- 
face are based on studies into the behavior of a single 
particle at the interface. The stability of an adsorbed 
colloid and the manner in which it attaches to the inter- 
face gives insight into the way particles act at higher con- 
centrations. The stability of colloids at an interface was 
already considered by Pieranski [JJ, who studied the ad- 
sorption free energy based on surface tension arguments. 
This ground breaking work was built upon to encompass 
effects, such as line tension [H [5], capillary rise [SI [7], 
surface deformation due to gravity [5], surface hetero- 
geneities [3 EH [UJ, and electrostatic effects O IT5] . 
The influence of particle shape on colloid adsorption has 
also been considered, for instance, ellipsoidal rods and 
platelets |H [B], and more complex shapes as well [T3] , 
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Nevertheless, there are still many unanswered questions 
concerning the adsorption of a single particle at an inter- 
face. 

To the best of our knowledge only one theoretical study 
has been undertaken into the effects of anisotropic parti- 
cles adsorbed to the interface as a function of the parti- 
cle's orientation [TS]. Most studies have been limited to 
several mathematically convenient particle orientations, 
namely parallel or perpendicular to the interfacial nor- 
mal SI [11 El [H HOI HH CH HZ]- These orientations are 
also found in experimental systems [SJ [HJ [THJ [THl HH1 HI] 
and therefore the current theoretical descriptions give in- 
sight into the behavior of the particles. However, these 
insights are constrained to particles that remain in one 
of these orientations. Therefore, these theories cannot be 
used to analyze the mechanisms by which colloids end up 
in these orientations, or why these particular orientations 
are preferred over other orientations. 

Studying the free energy associated with the adsorp- 
tion of an arbitrary shaped colloid with contact angle 
surface patterns is quite involved, especially when the 
colloid is allowed to have an arbitrary angle with the in- 
terface. We first examine homogeneous uni-axial convex 
colloids and formulate the adsorption free energy. Deter- 
mining this adsorption free energy proves to be techni- 
cally difficult for all but the most basic shapes. There- 
fore, we introduce a numerical technique, which we refer 
to as "triangular tessellation" , to evaluate the adsorption 
free energy. The accuracy of this technique is verified 
by comparison with semi-analytic results for ellipsoids, 
cylinders, and spherocylinders. These semi-analytic re- 
sults are derived by methods similar to those used in 
Ref. [IH]- We improve some of these results. Further- 
more, we extend the semi-analytic results of Ref. [TS] to 
a wider class of particles. Finally, we formulate a theo- 
retical description and present a numerical technique to 
handle non-convex colloids with surface patterns. 
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In conclusion, we introduce a new numerical scheme 
to determine the adsorption free energy of non-convex 
particles with or without surface patterns at the liquid- 
liquid interface, which has many advantages, as regards, 
applicability, stability, and generality, over semi-analytic 
techniques used so far. More detailed studies based on 
this new technique will be presented elsewhere |22j . 



II. METHOD 
A. Theoretical Considerations 

We consider a planar oil-water interface separating two 
homogeneous half spaces of oil and water, and a solid 
uni-axial convex colloid adsorbed at this interface. We 
focus here on an oil- water interface, but we note that any 
liquid-liquid interface can be considered, and to some ex- 
tent the theory is valid for liquid-gas interfaces as well. 
For simplicity, capillary effects due to the presence of the 
colloid at the interface are neglected. The coordinate 
frame is chosen such that the normal of the interface is 
along the z-axis. The position or depth of the interface 
with respect to the center of the particle, at which the ori- 
gin of the system is located, is denoted by depth z, which 
can be both positive and negative. For convenience we 
assume that the rotational symmetry axis of the particle 
is oriented in the xz-plane. The half space above the in- 
terface is called medium 1 (Mi) and the half space below 
the interface is called medium 2 (M2), see Fig. [l] The 
angle between the colloid's rotational symmetry axis and 
interfacial normal is denoted by <fi £ [0, ir/2]. Henceforth, 
<p is referred to as the polar angle. 

There arc four surface areas with corresponding surface 
tensions, which contribute to the adsorption free energy 
of the colloid: (i) the surface area of the colloid above 
the interface Si , (ii) the surface area of the colloid below 
the interface 5*2, (iii) the surface area excluded from the 
interface by the presence of the colloid Si 2, and (iv) the 
total surface area of the interface (without adsorption) 
A. There is also a contribution from the contact line, of 
length L, where the three phases meet, i.e., Mi, Mi, and 
the colloid. Writing S for the total surface area of the 
colloid, the following relations are obtained S — Si + S2 
and Si,S 2 G [0,S]. 

These surface areas and the corresponding surface ten- 
sions, together with the contact line length and corre- 
sponding line tension give rise to an adsorption free en- 
ergy. Such an adsorption free energy was first considered 
by Pieranski [1] and later extended to accommodate line 
tension, e.g., Refs. |H [T6]. For a specific configuration, 
characterized by z and <f>, this can be written as 

V(z, 4>) = 7l2 (A - S 12 ) + 7lc 5i + 72c S*2 + tL, (1) 

where 712 is the M1-M2 surface tension, 71 c is the M\- 
colloid surface tension, 7 2c is the M 2 -colloid surface ten- 
sion, and t is the line tension. Note that we have dropped 




FIG. 1: A colloid (ellipsoid) adsorbed to the interface lo- 
cated at depth — z, measured from the center of the colloid 
(0,0), dividing medium 1 (Mi) and medium 2 (M 2 ). The top 
frame shows an :rz-plane cross section of the colloid and in- 
terface. The polar angle between the interfacial normal and 
the rotational symmetry axis of the colloid is denoted by <f>. 
The interface has total area A with corresponding interfacial 
surface tension 712. The surface area of the colloid above 
the interface is denoted by Si with 7 i c the Mi-colloid surface 
tension and the surface area of the colloid below the inter- 
face is denoted by 5*2 with 72c the M2-colloid surface tension. 
The presence of the adsorbed colloid causes an area S12 to 
be excluded from the interface, indicated by the dashed curve 
in the xz-view. In the bottom frame this excluded surface 
area is shown from an xy-view, i.e., the region enclosed by 
the dashed line. This dashed curve is the contact line, with 
length L and corresponding line tension r. The outline of the 
colloid (solid curve) is also included. 

the z and </> dependence of S\(z, <f>), 6*2(2, 0), Si2(z, (f>), 
and L(z, <$) to lighten the notation. The tensions (surface 
and line) are system parameters. 

The systems to which Eq. can be applied, have neg- 
ligible capillary or gravitationally induced interfacial de- 
formation, i.e., they have a flat interface and insignificant 
flotation force. The validity of such an assumption for 
colloidal systems can be studied using the Bond-number 
parameter, which gives the ratio between gravitational 
and surface tension induced effects on an adsorbing par- 
ticle. Equation holds in the zero Bond-number limit 

Bo=^^0, (2) 

712 

where Bo is the Bond-number, g is the gravitational ac- 
celeration, Ap is the density difference between Mi and 
Mi, and R is the "characteristic" length scale associated 
with the particle. This result follows from conclusions 
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in Refs. [SJ [T7]. The length scale R is not entirely well- 
defined in the case of acicular (non-spherical) colloidal 
particles, but can be taken to be the mean radius of cur- 
vature. For colloids we estimate R to be at most 10~ 5 m. 
The surface tension 712 should be at least 1CP 3 N m _1 
for a liquid-liquid interface in the absence of surfactants, 
and the density difference can be at most of the order 
10 3 kg m~ 3 for physically reasonable systems. The Bond 
number is thus estimated to be Bo <C 1CP 3 . Hence it is 
safe to take the zero bond number limit [17], which im- 
plies insignificant interfacial deformation due to gravity. 

Capillary effects through immersion forces, anisotropy, 
and electrostatic interactions, cannot be so easily esti- 
mated. In Ref. [16 , capillarity is discussed in the context 
of inter colloid interaction on the interface, however, the 
effects on a single particle's adsorption are not included. 
The theoretical validity of neglecting capillary deforma- 
tion can be analyzed via de Young-Laplace equation, as 
done for a single orientation of ellipsoidal particles in 
Ref. [5J. This analysis and the consequences for compar- 
ison with experimental systems, pertaining to the theo- 
retical approach outlined in this paper, is left for future 
study 

It is customary to define the adsorption free energy 
with respect to a reference point. The shifted adsorption 
free energy F{z, <f>) is introduced by modifying V(z, <j>) in 
such a way that it is zero when the colloid is completely 
immersed in Mi, i.e., by subtracting 712A + ji c S. This 
yields 



F(z, 



(7xc - 72c)(<5i - S) - 712-512 + tL. (3) 



Often, the term "shifted" is ignored when referring to 
F(z,(j>). Only systems with Mi ^ M 2 are considered, 
hence 712 ^ 0. Therefore it is possible to write 712 cos 9 = 
7ic — 72c, where the contact angle 9 is introduced via 
Young's equation [23]. Using this definition Eq. ^ is 
rewritten to 



F{z,<j>) = 1i2 [{Si-S)cos9-S 12 ] + tL. 



(4) 



Note that the contact angle 9 is a quantity, which de- 
pends on the physical properties of the three components 
present at the interface, whereas the polar angle <f> is a 
degree of freedom. 

Dividing the adsorption free energy by 712s 1 (712 7^ 
0) and writing z = z* \J a 2 + 2b 2 , with a the rotational 
symmetry semi-axis, b the perpendicular semi-axis, and 
m = a/b the aspect ratio, the following dimensionless 
adsorption free energy is derived 

F< ^- = cose(ri-l)-ri 2 +T% (5) 
712J 

where r\ = Si/S and r%2 = S12/S are surface area ratios, 

r* = (6) 
712V0 

is the dimensionless line tension, and I = L/yS is a 

dimensionless contact line length. The value yja 2 + 2b 2 



is the length of the semi-diagonal of a rectangular beam 
with sides 2a x 2b x 2b. Two inequalities, namely < 
ri(z*,<f>) < 1 and < r 12 (z*,cf)) < 1, hold for any value 
of z* and 4> values. Note that f(z*, </>) implicitly depends 
on cos 9, m, and the shape of the particle. 

The dimensionless adsorption free energy, Eq. |5|, is 
scale invariant, i.e., it is independent of the size of the col- 
loid. For r* = our results for Eq. (|5| hold for any size 
colloid. For r* 7^ the results for Eq. (|5| can be trans- 
lated back to any colloid size, under the condition that for 
each size the line tension should be scaled. In this paper 
we focus on Eq. ^ rather than Eq. (4j), since our goal is 
to develop a general method to describe the adsorption of 
arbitrary colloid shapes at an interface. Therefore, semi- 
analytic and fully numerical results are compared for the 
dimensionless adsorption free energy. 

Before we present our numerical technique, we first 
introduce some quantities which prove to be useful in 
describing the results. The location of the adsorp- 
tion free energy minimum in Eq. ^ is denoted by 
(z* d , </> a d), which is referred to as the adsorption orien- 
tation. The corresponding adsorption free energy reads 
/ a d = f(z* d , 0ad)- Note that there may be multiple min- 
ima, in which case there can be meta-stable adsorption 
orientations. Several minima are labeled with a subscript 
i = 1,2, ... , where the deepest minimum has the lowest 
index. When there are two or more minima with equal 
adsorption free energy, we label them arbitrarily. 

For a given <j6, z*j et (<j6) is defined as the positive value 
of z* for which the interface just touches the top of the 



particle. The colloid is detached when z* < 



^det 



(4>) or 



z* > z dct (4>). The quantity z^ lin (0) is defined as the value 
for which f(z*,<j)) assumes its minimum as a function of 
z* for a given </>. The corresponding adsorption free en- 
ergy is denoted by / min (0) = /(CW> Note that 
it is a priori not excluded that for a given (f> the equi-0- 
curve has two or more (meta-stable) minima. This can 
correspond to multiple z^ in (0) curves in the free energy 
landscape, running "side-by-side" in the <j> direction. Of- 
ten we will write z det and z^ in for z* dct {<f>) and z^ in {4>) 
respectively, taking the polar angle dependence to be im- 
plicit . 

In summary, we have described our theoretical model 
for a colloidal particle at a planar interface based on free 
energy arguments in the zero Bond number limit. This 
model can be straightforwardly generalized to encompass 
more complex colloidal properties such as surface pat- 
terns, see Appendix A. Including the effects of gravity 
and capillary interfacial deformation is, however, sub- 
stantially more involved and therefore not pursued here. 



B. Numerical Approximation Scheme 

Determining the dependence of Si, S 2 , S12, and L on 
z* and <p is highly non-trivial in general, and deriving 
analytic expressions is unpractical, if not impossible, for 
all but the simplest particle shapes and orientations, see 
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Appendix B. To analyze colloids adsorbed at an inter- 
face the following numerical technique is employed. The 
surface of the colloid is bijectively parametrized by two 
angles, namely a± € [0,2tt] (azimuthal) and ev 2 € [0,7r] 
(polar) . A paramctrization can for instance take the form 

(r{a,\ , a.2) cos ol\ sin ev 2 \ 
r(ai, a 2 ) sin ax sina 2 , (7) 
r(ai, a 2 ) cosa 2 / 

where r(a.\,a.2) is some radial function, but many other 
forms are imaginable. The strip [0, 2tt] x [0, tt] is divided 
into triangles, the vertices of which are mapped onto the 
surface of the particle by means of the parametrization 
P(ai,a 2 ), see Fig. [2] A mapped triangle is formed be- 
tween the vertices of a corresponding triangle in the strip 
after P has acted on them. From now on, the object, on 
which the strip's triangle mesh is mapped, is referred to 
as being tessellated with triangles. The above method 
of modeling a 2D or 3D object by triangles (more gen- 
erally polygons) is well-known in computer science and 
has been successfully applied to various surface tension 
problems in physics [MJ [55] . 

The surface area of the colloid is now approximated 
by summing the surface areas of the mapped triangles. 
Suppose that the vertices of a mapped triangle are given 
by x, y, and z, then its surface area is given by a sim- 
ple cross-product |(z — x) x (y — x)|/2. This procedure 
can yield "in principle" arbitrary precision by sufficiently 
refining the triangular mesh. It should be noted that 
depending on the parametrization some triangles have a 
vanishingly small or zero contribution to the surface area. 
For example, in the case of a sphere several vertices co- 
incide resulting in degenerate triangles at the poles, see 
Fig. [2] Note that this mapping is not bijective, but only 
on a set of which the image has measure zero. 

The method described above can be amended in the 
following manner to enable calculation of Si, S 2 , S%% 
and L. Suppose that the tessellated object is inter- 
sected by a plane, then some of the triangles which 
compose the object lie above it and others below it. 
Let A 1 denote the set of triangles which lie strictly 
above, Ai the set of triangles which lie strictly below, 
and A p the set of triangles which intersect the plane or 
touch it. The surface of the colloid is approximated by 
S = Ev* A T,* + EyjA|_,J + J^yk A p,k, where the tilde 
indicates that this is an approximation, i, j, and k are 
indices, and the notation for an element in a set doubles 
as the notation for that triangles surface area. 

Each intersected triangle A p j is divided into three sub- 
triangles A rj i, A Sj i, and A t ^ in the manner indicated in 
Fig. [3j Two of these lie on one side of the plane, and one 
on the other. Applying this technique to all triangles 
in A p a set of partitioned triangles A p is obtained, of 
which the members only have some vertices in common 
with the plane and do not intersect it. Let A Pi j and A p ^ 
be the sets of triangles in A p which lie above and below 
the interface respectively, and let Aj = Aj U A p ^ and 
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FIG. 2: Example of a parametrization P (01,02) mapping 
the vertices of triangles in the strip (01,02) G [0, 2ir] x [0, tt] 
to points on a sphere. In between these points "mapped" 
triangles are formed corresponding to the original triangles on 
the strip. The points at the poles are degenerate, all vertices 
with 02 = and 02 = tt coincide at the corresponding pole. 
Only the front half of the sphere is tessellated for clarity. 




FIG. 3: Example of a triangle A Pi ; which intersects the in- 
terface (dashed line). Such a triangle can be cut into three 
pieces, A r ,i, A Sl i, and At.i as indicated above. In this case, 
the first piece lies above the interface and the second and third 
piece below it. The two new vertices lie on the interface. 



A 1 = A 1 U A Pj |. Using these sets of triangles the surface 
areas Si and S 2 are approximated in the following way 

Si = E A T,» (8) 

Vi 

£2 = £A W . (9) 

Vi 

Note that by virtue of this technique the equality S — 
Si + S 2 still holds, which can be used as a consistency 
check. 

From the set of triangles A p , the points where the 
plane intersects the original tessellation are extracted. 
These points form a two dimensional data set which ap- 
proximates the surface area that is cut out of the interface 
by the presence of the colloid, see Fig. [4j from which S12 
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FIG. 4: Graph (a) shows an example of the boundary points of 
the intersection of a colloidal cylinder with a planar interface, 
which are obtained using the triangular tessellation scheme, 
together with its convex hull (full curves). In (a) all points 
are shown, in (b) only those on the convex hull, with which 
the surface area S12 is computed through a trapezoidal inte- 
gration scheme as shown for the area above the x-axis. The 
convex hull also serves to compute the length of the contact 
line L. 



and L can be computed. 

Typically the boundary consists of several hundred 
grid points, depending on the size of the triangular mesh. 
The approximate surface area S12 is obtained by means of 
a trapezoidal integration scheme, which is applied to the 
points above and below the :r-axis after sorting them by 
increasing x-coordinate. Here we assume that the colloid 
is convex. The excluded area S12 is always a connected 
set when the colloid is convex. For non-convex colloids 
the area excluded from the interface can consist of two 
or more disjoined pieces, e.g., for a dumbbell. To avoid 
such difficulties we restrict ourselves to uni-axial convex 
colloids and refer the reader to Appendix A for a more 
general algorithm. It should, however, be noted that in 
the case of a dumbbell the surface areas are parts of cir- 
cles and spheres and therefore this shape can in principle 
be handled analytically 22J. 

The trapezoidal integration scheme suffers from insta- 
bilities due to small numerical uncertainties, which po- 
tentially interfere with the sorting algorithm, as illus- 
trated in Fig.[4ji. These problems can easily be overcome 
by considering the convex hull of the data set, Fig. [4j), 
thereby eliminating such "noise" . Considering the con- 
vex hull does come at the price of reducing the number 
of data points. However, for several hundred grid points 
this effect is negligible. The convex hull of the bound- 
ary points is also used to determine L, the approximate 
length of the contact line L. 

In summary, we have introduced an explicit tessella- 
tion scheme to compute the surface areas Si, S2, and 



Oblate (a < b) 





Prolate (a > b) 






FIG. 5: Impression of the various colloidal shapes considered 
in this paper. The left column represents oblate particles 
(a < b, m < 1), the right column prolate particles (a > b,m > 
1). Note the difference in shape between oblate and prolate 
spherocylinders . 



Si 2 as well as the contact line length L explicitly, for any 
convex (uni-axial) colloid particle. This scheme can, how- 
ever, be straightforwardly generalized to more complex 
colloidal shapes with surface heterogeneities. Non-convex 
shapes and particles with surface patterns are discussed 
Appendix A. 



III. RESULTS AND DISCUSSION 

In this section we discuss the adsorption free energy 
landscape for several particle shapes, namely ellipsoids, 
cylinders, and spherocylinders, as shown in Fig. [5] We fo- 
cus on these three types of particle, as they are frequently 
used to model colloidal platelets and rod-like colloids in 
theoretical work and computer simulations. In addition, 
these particle shapes can be either prolate (m > 1) or 
oblate (m < 1), and have relatively simple parameteri- 
zations. Note that in our model the spherocylinder has 
length a and width b. Contrary to traditional notation 
a includes the spherical end-caps for a prolate sphero- 
cylinder, while b includes the rounded side for an oblate 
spherocylinder . 

Only the dimensionless adsorption free energy, Eq. (|5| , 
is considered and the investigation here limits itself to 
two aspect ratios m = 1/4 and m — 4, one contact an- 
gle cos 9 = —1/2, and several values of r*. Negative 
values for cos 9 are used, since then the particle prefers 
Mi. However, there is no real difference between cos# 
and — cos 9, because it only implies interchanging Mi 
and M 2 in our model. That is to say, when (z* d , (j)^) is 
the location of a minimum for cos 9 > 0, then for — cos 9 
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there is a minimum at (— z* d , c/) a( j) if r* is fixed. It should 
be noted that the line tension can be both positive and 
negative [37] and can assume values in a range spanning 
several decades [55] HHJ 130] depending on the precise de- 
tails of the system. The line tension r* is therefore chosen 
arbitrarily, without reference to a specific experimental 
system. We have limited ourselves to several interesting 
configurations, to prove the accuracy of our method. 

An equidistant mesh of 200 x 200 to 400 x 400 vertices 
is employed, which via the various parameterizations is 
mapped onto a heterogeneous triangular tessellation. For 
future reference we will denote a 200 x 200-vertex grid 
as a 200 2 -vertex grid, for instance. This yields relative 
fractional accuracies in the range 10 -3 to 10~ 5 in S, Si, 
S2, 5i2, L, r\, ri2 and I, depending on the object pa- 
rameterized. These uncertainties are established using 
semi-analytic values of the surface areas and contact line 
length, given in Appendix B. The semi-analytic nature 
refers to the fact that one dimensional integrals need 
to be evaluated in order to obtain a value. More than 
5000 non-equidistant grid points are used in these calcu- 
lations to ensure a relative fractional uncertainty lower 
than 10 -6 (the magnitude is determined using grid re- 
duction). These semi-analytic results are independent of 
mesh size and triangular tessellation and can be used to 
test our method. For higher accuracy, meshes of 1000 2 
vertices are employed, although in most cases these re- 
sults are indistinguishable from the 400 2 -vertex mesh re- 
sults, or the semi-analytic results. We are therefore con- 
fident that the numerical scheme is sufficiently stable and 
can be applied to shapes for which we have not performed 
analytic verification. 



A. Ellipsoids 

We calculate the adsorption free energy f(z*,(f>) of 
an ellipsoid with aspect ratio m = 4 and contact an- 
gle cos (9 = —1/2. We use a 400 2 -vertex mesh for the 
tessellation and a 5000-point equidistant trapezium-rule 
grid for the semi-analytic approach. Figure [6^l shows 
f(z*,<fi) for t* — in a 3D representation as a function 
of z* and (p. In Figs. |6^b-f) we plot (^-sections of the 
free energy landscape for varying r*. The numerical re- 
sults established using the tessellation scheme presented 
above and the semi-analytic results obtained using the 
equations and techniques described in Appendix B agree 
within the line width of the curves, i.e., the relative un- 
certainty is lower than 10 -3 for all grid points. Graphs 
similar to those in Fig.[6]can be made for ellipsoids, cylin- 
ders, and spherocylinders of any aspect ratio, with any 
contact angle, and line tension. We have verified that the 
semi-analytic scheme in Appendix B and our tessellation 
scheme yield the same results. 

In the specific case of an ellipsoidal colloid with m = 4 
and cos6> = —1/2, see Fig. [6] we find that for negative 
values of r* there is a single minimum in all 0-sections 
of the free energy landscape. From these figures it can 



also be derived that there is in fact a single minimum 
in the adsorption free energy for r* < 0. That is to 
say, a single minimum at <f> — tt/2 and no meta-stable 
secondary minima. We will come back to this shortly. 
From Figs. [6](b-f ) it is clear that the single minimum in 
the ^-sections with — z det < z* < z dct for r* < vanishes 
with increasing r* . The minimum free energy for a given 
<j), is then / m in(0) = and this minimum is assumed 
when z* < — z det . In this case, the particle prefers to be 
detached from the interface and can move freely in Mi, 
where its adsorption free energy is lowest. 

For certain r* a 0-section can have two minima, for 
example r* = 0.05 and cj) = tt/2, see Fig. [6J 7 . Here there 
is an absolute minimum given by the detached state with 
/min( 7r /2) = and a meta-stable (local) minimum with 
z* 6 [— z dct , z det ], i.e., when the particle is adsorbed at 
the interface. For sufficiently positive r* the presence 
of this local meta-stable minimum in the ^-sections is 
dependent on the value of <j). There is a local minimum 
with z* G [— z dct , z det ] for 4> = tt/2 and 4> — 37r/8 when 
r* = 0.025 for instance, but this minimum is not present 
for the < 7r/4 sections, see Figs. |6jb-f). Conversely the 
minimum with z^ in G [— z dct , z^ c Jcan be the absolute 
minimum and the detached state (z* < —z dct ) a meta- 
stable minimum, see for instance the (^-section in Fig. [Of 
with t* = 0.025. 

The appearance of a local meta stable minimum with 
z* G [— z dct , z dct ] is in part related to the formation of 
"adsorption barriers" in the free energy. That is to say, 
when r* is sufficiently positive, the colloid has to cross a 
barrier to attach at the interface from an immersed state 
in either medium. From Figs. (6^b-f) we observe that 
positive values for the line tension r* give rise to these 
adsorption barriers. This behavior is most clearly visible 
in the r* = 0.025 and r* = 0.05 ((f) = 7r/2)-sections in 
Fig.[6f. The height of the barriers varies with the value of 
the polar angle and they become more pronounced with 
increasing r*. These barriers are quite intriguing, since 
they can prevent a particle form reaching its lowest free 
energy state, when it is initially introduced in its least 
favored medium. 

The above results agree with the findings in Ref. |15j . 
However, there are several differences between our results 
and those of Ref. [15] as the expressions in Appendix B 
and in Ref. [TS] do not agree completely. Apart from 
minor typographical errors, there is a problem with the 
definition of subdomains on which the equations hold, as 
well as the way in which quantities are made dimension- 
less. Because of notational differences these problems are 
not immediately obvious, but when comparing results it 
is clear that the adsorption free energy barriers induced 
by the line tension are far less pronounced in our case. We 
believe that the way in which t is made dimensionlcss in 
Ref. [T5] violates scale invariance, but from the descrip- 
tion given this is impossible to determine. Despite the 
discrepancies with previously established results, we are 
confident that our results are correct, since we have used 
two independent methods, which yield the same results 
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FIG. 6: The adsorption free energy landscape for an ellipsoid with aspect ratio m = 4, and contact angle cos# = —1/2. Graph 
(a) shows f(z*, 4>) in a 3D representation as a function of z* and <j> for r* = 0. Graphs (b) through (f) show ^-sections of this 
free energy landscape for = 0, tt/8, tt/4, 3tt/8, and tt/2 respectively, and for r* = -0.1, -0.05, -0.025, 0, 0.025, 0.05, 0.1. 
The central dotted line in each graph corresponds to a section of graph (a), i.e., r* = 0. The adsorption free energy (^-sections 
are shown for —1 < z* < 1, but can be extended with constant value for z < —1 and 0.5 for z > 1. Note the appearance of 
free energy barriers for r* > which must be crossed if the particle adsorbs to the interface. When r* > 0.05 the adsorption 
free energy does not have a minimum. In this figure numerical and analytic results agree within the line width of the curves. 
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sorption free energy f m i n {<t>) = /(^mim^) f° r ellipsoids 
with m = 4, m — 1/4, cos# = —1/2, and several 
r* corresponding to the choices in Fig. [6] The mini- 
mum adsorption free energy is / m in(</ ) ) under the con- 
straint that z* nin e [~ z det' ^detl- The trivial solution 
/min(0) = with z* < — z*j ct , when r* 3> 0, is repre- 
sented by a thick horizontal line. Note that this solution 
is not in the [— z^ et , z^ et ] domain. We have also indicated 
any meta-stable part of the adsorption free energy, i.e., 
/min(0) > 0. These minima are metastable with respect 
to the / m in('/ > ) = solution, for which z* < — z*j ct . This 
local minimum need not exist for all values of <fi, which 
results in the minimum adsorption free energy curves 
terminating when r* > 0.025 (to = 4) and r* > 0.0 
(to = 1/4) respectively, see Fig. Uj This is in agreement 
with the behavior of the ^-sections given Fig. [6jd-f). For 
r* < 0.025 all minimum adsorption free energy curves 
are stable. 

If there is a minimum adsorption free energy curve for 
a given r*, then we find that the (local) minimum of 
the free energy is located somewhere on this curve, per 
definition. For to = 4 the location of the the adsorp- 
tion minimum of / corresponds to </> a d = tt/2. Similarly, 
when rn = 1/4 the adsorption free energy / is minimal 
for a d = 0. This observation proves to hold in general 
for ellipsoids, because the qf> a d is determined solely by 
the aspect ratio [25]. We alluded to this result earlier, 
in describing Fig. [6] but the presence of a single mini- 
mum and the monotonicity of the minimum adsorption 
free energy curves is much more evident in Fig. [7] To 
summarize, the effect of positive line tension is found to 
be two fold. Firstly, it gives rise to an adsorption barrier, 
which must be overcome for particles to attach to the in- 
terface. Secondly, it destabilizes the adsorption of a col- 
loid, by reducing the depth of the free energy minimum 
at adsorption, even to the point that it is either meta- 
stable or non-existent. The angular dependence of this 
stability reduction can be attributed to the anisotropy in 
the particle shape. 



* (it) 



FIG. 7: The minimum adsorption free energy f m in(4>) f° r 
ellipsoids with aspect ratio m = 4 (a), m = 1/4 (b), 
cos# = —1/2, and several values of r*. The thick line in- 
dicates / m in(0) = 0, which is the minimum free energy when 
r* 3> 0. Note that the meta stable part of the adsorption 
curves is shown as well, i.e., /min(</>) > 0. We only concern 
ourselves with the presence of a minimum for a certain <j) 
when z* € [— z^ et , z^ et ], not with its stability. Depending on 
the value of r*, a minimum in the 0-sections need not be 
present along the entire 0-range, also see text. For r* =0.1 
in graph (a) there is no minimum adsorption free energy, in 
correspondence to the results from Fig. [6] Again the agree- 
ment between numerical and semi-analytical results is within 
the line width of the curves. 

In addition to giving ^-sections, we have determined 
for which z^ in the adsorption free energy f(z*, 4>) is min- 
imal for a given <fr. Figure [7] shows the minimum ad- 



B. Cylinders 

We only study the minimum adsorption free energy 
as a function of the polar angle <fr for cylindrical par- 
ticles and several t* when cos6> = —1/2. Figure [8] 
shows the corresponding curves / m in(</ ) ) for cylinders with 
to = 4 (a) and to = 1/4 (b). Again we find that for 
t* > 0.1 there is only one minimum in the free energy 
with z* lin < — z^ ct , namely when the colloid is completely 
immersed in Mi . It can be shown that this trend holds in 
general. When comparing Fig. [8^, with the results given 
in Ref. |15j . there is no correspondence. This lack in ac- 
cordance can be attributed to the fact that the equations 
presented in Ref. [15] do not sample all possible orienta- 
tions of the cylinder with respect to the interface. Again 
our results have been verified by comparison with the 
semi-analytic result. The agreement is better than that 
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FIG. 8: The minimum adsorption free energy curves / m in(0) 
with cos# = —1/2 and several values of r* for a cylinder 
with m — 4 (a), and m = 1/4 (b). The thick line indicates 
/mm(</>) = 0, which is the stable minimum when r* 0. The 
agreement between numerical and semi-analytical results is 
within the line width of the curves. Note the appearance of 
a secondary minimum in /mm (</>)• The kink in the minimum 
free energy curves (graph b) for the oblate cylinder is caused 
by the sharp edge of the cylinder, see text. 



the line width of the curves in Fig. [8] can show. 

The cylindrical particles have two noticeable differ- 
ences with their ellipsoidal counterparts. Firstly, we find 
that two minima can appear in a single minimum ad- 
sorption free energy curve. Whereas for ellipsoids these 
curves have monotonic properties, and the free energy 
therefore has only one minimum, the cylindrical curves 
can have two minima. For the configurations considered 
the absolute minimum is located at <fi&d = tt/2 for to = 4, 



i.e., the particle lies fiat on the interface. The meta-stable 
minimum is at <f> a< \ = 0, i.e., the particle is perpendicu- 
lar to the interface. The latter corresponds to a cylinder 
which only has one of its disk shaped end-caps flush with 
the interface and the rest of its surface in M\. The free 
energy gained by excluding a disk from the interface and 
not having any other part in contact with M 2 , which is 
energetically unfavorable, can be sufficient to generate a 
local meta-stable minimum. It is also possible to choose 
parameters such that this "odd" configuration has the 
lowest adsorption free energy and is therefore stable. Sec- 
ondly, we also observe the presence of a kink in some of 
the minimum adsorption free energy curves, see Fig. [8J). 
This feature is related to the sharp corners of the cylin- 
der itself, which gives rise to ridges in the adsorption free 
energy landscape. For the two configurations considered 
here, only the minimum free energy curves for to = 1/4 
follow parts of these ridges, and consequently a kink ap- 
pears. For the prolate cylinder the minimum is found 
away from the ridges in the free energy landscape and 
the kinks do not appear in the minimum adsorption free 
energy curves. 



C. Spherocylinders 

For spherocylinders we again concentrate on determin- 
ing for which z^ in the adsorption free energy f(z*,4>) 
is minimal for a given <f>, under the constraint that 
z mhi ^ [— ^dct , ^dotl- The minimum adsorption free en- 
ergy curves f m in(4>) for spherocylindrical particles with 
to = 4 (a), m — 1/4 (b), cos = —1/2 are given in 
Fig. [9] for several values of r*. Note that for sphero- 
cylinders these curves are similar to those for ellipsoids. 
The agreement between numerical and semi-analytic re- 
sults is again better than can be appreciated from the 
line width. For the semi-analytic results a grid of 8000 
non-equidistant points was used to evaluate the one di- 
mensional integrals using Aitken's method |31j . also see 
Appendix B. This yields a relative error of 10~ 7 based 
on grid reduction. The numerical tessellation scheme is 
based on 1000 2 vertices to obtain a relative accuracy of 
10~ J per point or better, when compared to the semi- 
analytic curves. 

There is only one minimum per minimum adsorption 
free energy curve, as the particle shapes are smooth and 
in that sense more similar to ellipsoids than to cylin- 
ders. For spherocylinders, (j) at \ is also completely deter- 
mined by the aspect ratio m, i.e., a( j = for m < 1 
and a d = vr/2 for to > 1 [22 j. This property further 
distinguishes smooth convex uni-axial particles, e.g., el- 
lipsoids and spherocylinders, from non-smooth particles, 
e.g., cylinders. However, the mechanism behind adsorp- 
tion to the interface is quite subtle, depending not only 
on to, S, cos 8, and t* , but to a large degree on the 
shape and surface patterning of the colloid. A more de- 
tailed study on the effect of particle shape and surface 
patterning will be the topic of future work |22j and will 
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FIG. 9: The minimum adsorption free energy curves / m in(0) 
with cos# = — 1/2 and several values of r* for a spherocylin- 
der with m = 4 (a), and m = 1/4 (b). The notations are 
analogous to those in Figs. [7] and [5] Again the agreement be- 
tween the two methods used to generate the results is better 
than the line width of the curves. Note that in both graphs 
there is only one <^> a d per minimum curve. 



be presented elsewhere. 



IV. CONCLUDING REMARKS 

We have presented a numerical framework to deter- 
mine the adsorption free energy, in terms of surface and 
line tension contributions, of a non-convex patterned col- 
loid adsorbed at a flat interface. This framework is a 
natural extension of established theoretical models, e.g., 
Refs. [TJ HI US]. Because of the complexities which 



arise when determining the adsorption free energy for 
complex colloidal shapes, a numerical technique based 
on triangular tessellation is developed. The accuracy 
and stability of this method has been extensively veri- 
fied for convex uni-axial colloids, such as ellipsoids, cylin- 
ders, and spherocylinders. This analysis was performed 
by comparing semi-analytic results with the results pro- 
duced by our tessellation technique. The expressions 
used to obtain the semi-analytic results are given in Ap- 
pendix B, and thereby amend previous results |15J . Ex- 
act correspondence is found between the numerical and 
semi-analytical results, and the established analytic val- 
ues for specific configurations of colloids adsorbed at a 
flat interface [4]. 

Despite the fact that all of the presented results can 
be derived semi-analytically, the numerical method pre- 
sented here has substantial merits. On the one hand, 
the semi-analytic results can be calculated numerically 
substantially more quickly than the triangular tessella- 
tion results. Only the value of one dimensional integrals 
needs to be approximated for these. On the other hand, 
the process of determining analytic expressions can be 
very labor intensive and may result in equations which 
are incorrect. Therefore, the numerical scheme can be 
employed to verify the equations derived by analytic 
means. However, verification is not the true strength 
of the numerical technique. Determining an appropriate 
parametrization and corresponding tessellation is much 
less involved than obtaining the analytic expressions, and 
hence the numerical method is more suited to quickly de- 
termine the adsorption free energy landscape of a range 
of different shapes and surface patterns. 

It should also be pointed out that, once implemented, 
the triangular tessellation scheme will have the same nu- 
merical limitations for any system it is applied to. By 
identifying these limitations, they can be avoided, mak- 
ing the method very robust. Semi-analytic results are 
different in this respect, since they suffer from a lack of 
uniformity in the numerical techniques that are needed. 
The diversity of ID and 2D integrals, which require nu- 
merical evaluation, obtained by studying general systems 
is limitless. The numerical evaluation of these should be 
scrutinized on a case by case level, a numerical integra- 
tor that works for one integral, is not necessarily suited 
for another. In our experience, the merits of the tri- 
angular tessellations scheme, greatly outweigh those of 
the semi-analytic approach for generality, stability, and 
applicability. Hence, our method based on such a tri- 
angular tessellation scheme, can be used to examine the 
physics behind interfacial adsorption of general colloidal 
particles as well as to make a link between theory and 
experiments |22j . 
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APPENDIX A: EXTENSION OF THEORY AND 
NUMERICAL SCHEME TO NON-CONVEX 
PATTERNED PARTICLES. 

1. Theoretical Model 

In this section we extend our theory to non-convex col- 
loids with surface patterning. We consider a colloid, of 
which the surface is divided into several areas with dif- 
ferent liquid-solid surface tension properties, as is illus- 
trated in Fig. [TO] The surface patterning, other than the 
most basic, breaks the rotational symmetry properties of 
our system. Hence, an arbitrary non-convex patterned 
particle is described by an angle to in addition to the two 
parameters z and <j> used earlier, see Fig. [10] Suppose 
an arbitrary axis is fixed through the center-of-mass of 
the particle, at which the origin of the system is located, 
then 4> is the angle between this axis and the interfacial 
normal. The angle w describes rotations around this axis. 
Rotations by u are well defined, if they are with respect 
to some initial orientation. This initial orientation can, 
however, be arbitrarily chosen. 

Let the particle be partitioned into n patches Pi, with 
the index i — 0,1, ... ,n, see Fig. [10] Each patch has a 
specific patch-medium surface tension 7i c . (Pj-Mi) and 
72e; (P1-M2), in addition to a line tension Tj. Let Si i 
be the surface area of patch Pi in Mi, S 2i the surface 
area in M 2 , and L t the length of the contact line. These 
three quantities can be any value between and their 
maximum, depending on the orientation of the colloid. 
The adsorption free energy of the colloid is given by 

n 

V(z,<j>,uj) = Jl2{A-Si2)+y~] (TIcA + l2c z S 2 , + nLi) . 

i=0 

(Al) 

We define Si = S^ + S 2i , such that the total surface area 
is S = Yl7=o As before we can set the adsorption free 
energy to be zero in Ml, by subtracting 



7i2^+y^7ic,^ 



(A2) 



from V(z, 4>,uj), to obtain the shifted adsorption free en- 
ergy 

n 

F(z, (f>,uj) = y~] ((71c. - 7 2ci )('5i j - Si) + TiLi)-^ 12 S 12 . 

z=0 

The following quantities are introduced analogous to the 
theoretical description given earlier, 712008^ = 7 lc . — 
72ci via Young's equation, = Si/S, r\ f = S\JS, 
Tvi = Si 2 /S, k = Li/s/S, and r* = n/j^y/S. Using 




FIG. 10: A non-convex surface-patterned colloid adsorbed to 
an interface in schematic representation (top) and the same 
colloid's surface patterning (bottom). The orientation of the 
colloid is described by three quantities: (i) the depth z; (ii) the 
angle <f), which some arbitrary but fixed axis through the col- 
loid's center-of-mass makes with the interfacial normal; (iii) 
the angle uj which gives rotations around the chosen axis. The 
angle uo is defined with respect to a predetermined initial con- 
figuration. The origin (0, 0) is located at the center-of-mass 
of the colloid. Again A is the total (macroscopic) surface area 
of the interface, S12 the area excluded from the interface by 
the presence of the particle, and 712 the medium 1 (Mi) - 
medium 2 (M2) surface tension. The individual patches are 
labeled by Pi (bottom), where i = 0, . . . n. Each patch has a 
Pi — Mi surface tension 71^ , a P, - M2 surface tension 72^ , 
and a line tension r^. The surface area of Pi in Mi is given 
by Si t , the surface area in M2 by 5^, and the contact line 
length by Li. 



the above definitions Eq. (A3 1 is reduced to the elegant 
dimensionless from 

n 

f(z*,(f>, u) = (cos 9 l (r u - n) + r*k) - n 2> (A4) 

i=0 

where the shifted adsorption free energy has been divided 
by 712 S. Here z is made dimensionless by introducing 
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z* = z/TZ, with TZ the radius of the smallest sphere which 
encloses the particle. 

Note that in our model we have glanced over one de- 
tail, namely that it is possible to have a four-phase con- 
tact line on the boundary of two patches, if the interface 
coincides (partially) with this boundary. We define the 
line tension associated to this four-phase contact line to 
be the average of the line tensions of the two patches, 
for mathematical convenience. The contribution to the 
adsorption free energy is the length of this four-phase 
contact line times the averaged line tension. To the au- 
thors' knowledge, little is known about the properties of 
such four-phase line tensions. There have been studies 
into four-phase contact lines but the line tension is 
not included. We believe that averaging is not unrea- 
sonable. However, multi-phase line tensions and point 
tension contributions to the adsorption free energy cer- 
tainly merit further investigation. Sharp features in the 
colloid, such as cusps and facets, and their possible ad- 
sorption free energy contributions due to stresses induced 
on the interface by surface exclusion are also not consid- 
ered here. 



2. Improved Numerical Scheme 

The numerical scheme, based on triangular tessellation 
as described in the main part of this paper, is suited to 
handle non-convex shapes, with the exception of the sur- 
face area excluded from the interface. In this paragraph 
we extend the numerical scheme to handle non-convex 
patterned particles. Note that it is not excluded that 
these particles have handles. 

First the object is tessellated with triangles using a 
suitably chosen parametrization. This tessellation obeys 
the following rules, (i) Patch boundaries are approxi- 
mated by triangle edges, consequently, a single triangle 
has a single set of surface properties, (ii) Each triangle is 
labeled according to the patch it is in, with the label Pi. 
(iii) The direction of the surface normal of each triangle 
is known, and is required to point outwards from the par- 
ticle, (iv) A sufficiently large number of small triangles 
is used where the surface of particle changes abruptly, 
cither via a large gradient or a cusp-like feature. 

Determining the approximated total surface area S is 
analogous to the procedure outlined before. Similarly, 
the surface area of patch i, Si, is approximated by sum- 
ming the surface area of triangles with label Pi , yielding 
Si. The surface area of a patch above the interface S± i is 
obtained by partitioning the triangle mesh in the man- 
ner described earlier. The approximated surface area is 
denoted Si 4 . The area can also be determined in 
this way. The equality Si = Si i + also holds in ap- 
proximated form Si = Si i + S2i and can be used as a 
consistency check. 

The calculation of S12 and Li is, however, a little more 
involved. In partitioning the triangles, two points are 
obtained for each triangle, if the sides of the original tri- 



angle are intersected by the interface. These two points 
span a line segment, which is oriented via the normal of 
the triangle. That is to say, after partitioning a set of line 
segments, say A, is obtained, see Fig. [TTfc . The members 
of A, say Aj with i an index, are encoded with infor- 
mation on the location of the particle. This encoding is 
as follows. For a triangle intersected by the interface at 
an angle, the triangle's normal is projected onto the in- 
terface and normalized. This unit vector is referred to 
as the directional (unit) vector, because it gives orienta- 
tion to the line segment. A triangle which lies flush with 
the interface needs to be special cased, here the word 
flush indicates that all vertices are located in the plane 
of the interface. The directional unit vectors for each of 
its sides point outward, the implementation of which is 
trivial. Each vector Aj has 7 components, two give the 
xy-location of the starting point of the line segment, two 
the end point of the segment, two give the direction of the 
unit vector, and one gives the patch it is associated to. 
All directional unit vectors obtained in this way, point 
outward from the colloid (at least locally). 

To determine §12 the set of line segments A needs to 
undergo several refinement steps first. Consider all in- 
stances of a A,-, for which there is a Aj, which has the 
same line segment coordinates, but not necessarily the 
same directional vector, e.g., see Fig.fTTfc. The situations 
in which there are two "overlapping" segments are the 
following, (i) When one of the sides of a triangle is flush 
with the interface (all points off that side are located in 
the plane), this gives one line segment. The second seg- 
ment is given by the triangle which shares that particular 
side with the original triangle. Both of these segments 
have the same directional vector, (ii) When an entire tri- 
angle is flush with the interface all three sides contribute 
a line segment. These segments need not necessarily be 
a part of the boundary of Si 2, or equivalently not a part 
of the contact line, they may be internal, see Fig. [TT^ .. If 
one of the sides of the original triangle is an interior side, 
the adjacent triangle, must also be entirely flush with 
the interface. Otherwise, it would not be an internal seg- 
ment, i.e., lie in the interior of Si 2. This ajacent triangle 
gives a second instance of the line segment, for which the 
directional vector is opposite to that of the first segment, 
sec Fig. 11 1. To eliminate unnecessary segments Aj is 
removed from A if Aj and A^ have the same directional 
vector. If, however, the directional vectors have opposite 
sign both instances are removed, since then these are in- 
terior segments. By subjecting each element in A to this 
procedure a new set A is formed. This set contains only 
segments which are a part of the boundaries of Si 2, sec 
Fig. [lip. 

The set A is subdivided into closed loops. A loop is 
determined by choosing a segment in A and adding its 
neighbors recursively, until no more new neighbors can 
be added. This procedure is illustrated as follows. Let 
Aj be the starting segment. Then after one iteration 
we obtain the sequence A,_i — Aj — A J+ i, and after m 



iterations A;_ r 



-Aj_i— Aj— A, 



The 
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FIG. 11: Illustration of the loop reduction procedure for 
the approximated interfacial cut-out 5i2 corresponding to a 
non-convex colloid. We chose not to complicate the figure by 
composing S12 out of disjoined pieces, but the procedure for 
this situation is analogous. Box (a) shows the line segments 
obtained after partitioning triangles, i.e., the set A. Note that 
some sides have multiple directional unit vectors, indicated 
by arrows, associated to them. These are double instances of 
a line segment, which occur whenever triangles have one or 
more sides flush with the interface as is explained in the text. 
Box (b) shows A, the set of segments after removing all double 
and internal line segments. Box (c) shows the subdivision of 
A into three closed oriented loops. 



last neighbors to be added are either equal, i.e., Aj_ m — 
Ai +m , in which case only one is added, or have a common 
vertex, in which case the loop is also closed. We thus 
obtain a loop which is ordered by construction. This 
procedure is repeated until A is subdivided into loops. It 
is a priori not excluded in the above that a loop crosses 
itself, e.g., a lemniscate like structure. Crossover points 
are however easily located, by the fact that such a vertex 
will have attached to it an even number of line segments 
greater than two. All loops are subsequently cut into 
closed pieces which do not cross themselves, see Fig.pTp. 
Let these loops be denoted by Ci with i and index, then 
A = UviCi and n Vi Ci = 0. 

It is necessary to determine the type of closed loop, i.e., 
whether it is outward or inward. By construction loops 
cannot contain both outward and inward segments. For 
a loop where all the directional vectors point inwards the 
following holds. Any half-line, starting either end-point 
of a line segment in the loop, in the direction indicated 
by the directional vector, will intersect another line seg- 



ment in the loop. For a loop where all directional vec- 
tors point outward, there is at least one line segment for 
which one of the half-lines drawn through its end-points 
as before will not intersect another line segment in the 
loop. Numerically checking this criterion efficiently is 
not trivial, however, the maximum length over which a 
half-line needs to be checked for intersection is 27?., with 
TZ as in the definition of z* . After each loop has been 
labeled either "outward" or "inward" , the area enclosed 
by each Ci in A is easily calculated using a polygonal 
version of Green's theorem. In practice, it is seldom re- 
quired to use such a complicated scheme to determine the 
orientation of the loop. Any knowledge on the possible 
interfacial cut-outs can be used to make the algorithm 
more efficient. A dumbbell, for instance, only has out- 
wardly oriented loops, which can be easily derived from 
its symmetry properties. 

The line segments in a loop Ci define a set of 2D points 
in the plane, which are ordered by the ordering of the 
loop. Map these points onto three dimensional vectors 
a.j, where the first two components are xy-coordinates, 
the last component is zero, and j is an index. Let the 
set of these vectors be ordered according to the ordering 
imposed by Ci. If there are h — 1 distinct points which 
define the loop, let j = 0, . . . , n with = ao- The area 
of Ci is then given by 



3=0 



(A5) 



where the x -symbol indicates the cross product and the 
vertical bars the norm of the vector obtained by summa- 
tion. Note that this is indeed an adaptation of Green's 
integral theorem to polygonal shapes. It can be shown 
that Eq. ( |A5[ ) is only valid when there are no self inter- 
sections, which is why these needed to be eliminated first. 
Define S(Ci) to be +1 when the orientation of the closed 
loop is outward and —1 if the orientation is inward. The 
approximated surface area excluded from the interface 
by the presence of the particle is given by 



S 12 = Y,S(C l )A(C t ). 



(A6) 



Vi 



The calculation of the approximated contact line lengths 
Li is also possible from the set A, by summing over the 
lengths of line segments which have the same label Pi. 
Note that one needs to special case the instances when 
the contact line coincides with patch boundaries, as de- 
scribed in the previous section. 

In summary, we have presented a method that can be 
used to determine the free energy of an arbitrary pat- 
terned colloid adsorbed to a flat interface. This numerical 
scheme will be applied in future work [22 on for instance 
colloidal dumbbells and particles which have Janus-like 
patterning. There are, however, still open problems, such 
as point tensions, cusp-effects, and four-phase line ten- 
sions which merit further investigation. Finally, it should 
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be noted that in specific cases this scheme can be greatly 
reduced if properties of the possible interfacial cut-out 
shapes are known. 



APPENDIX B: ANALYTIC EXPRESSIONS FOR 
THE ADSORPTION FREE ENERGY 

In this Appendix we reproduce the analytic expressions 
for ellipsoids, cylinders, and spherocylinders used to ver- 
ify the accuracy of the triangular tessellation method. 
To keep Appendix B concise, derivation of the results is 
not included. The methods used are however analogous 
to those applied in [TS]- Only the expressions for S, Si, 
S12, and L are considered here. We aim to keep the for- 
mulation as general as possible, while at the same time 
showing the similarities and dissimilarities between the 
various shapes. 

Some expressions can be reduced in specific cases, i.e., 
several integrals for ellipsoids can be evaluated to give 
closed expressions in terms of standard functions. How- 
ever, such a reduction may result in the expression only 
holding for oblate and not for prolate particles or vice 
versa. Some integrals can be reduced using symmetry 
properties, this will not be done here in order to em- 
phasize similarities. In the case of a spherocylinder the 
strong difference in shape between oblate and prolate ne- 
cessitates differentiation between the two aspect ratios. 
A general equation which describes both types of sphe- 
rocylinder cannot be given. 

It is always implied that the following symmetry prop- 
erties 



S 2 (z,(f>) 

Si 2 (-z,4>) 
L(-z,<p) 



S — Si(z, i 
S — S±{z, i 
Si2(z,4>); 
L(z,4>), 



(Bl) 
(B2) 
(B3) 
(B4) 



are used to describe the system and speed up numerical 
calculation. The equations and variables considered are 
therefore only given on the {z, 0)-domain 



D = [0,oo) x [0, tt/2] • 



(B5) 



It will prove necessary to subdivide this domain into "dis- 
joined" pieces on which equations are defined. It can 
be shown that the equations defined on these subdo- 
mains change into each other continuously on the com- 
mon edges. It can also be shown that these equations re- 
duce to previously established results, Refs. [HQS], when 
= and (f) = 7r/2. We do not include these calculations 
here in the interest of briefness. 

Because of the large number of symbols required to 
formulate the expressions for S, Si, Si 2, and L, we are 
forced to recycle notations on a paragraph by paragraph 
basis. However, an attempt is made to use the same sym- 
bols for similar quantities as much as possible. The same 
holds for the definitions of the subdomains for the vari- 
ous species of particle. To further reduce the notation the 



dependence of variables on z and <f> is often implicit. In 
the case that a parameter has a different value on several 
subdomains, it is implied that any function depending on 
this parameter should be evaluated with the appropriate 
value. 

All particles considered analytically require the numer- 
ical evaluation of one dimensional integrals in order to 
calculate the various surface areas or the contact line 
length. Therefore, we refer to this method as semi- 
analytic. For ellipsoidal, cylindrical, and prolate sphe- 
rocylindrical particles a simple equidistant trapezoidal 
scheme can be implemented with relatively small numer- 
ical error. However, a more stable technique is required 
for oblate spherocylinders, due to divergences in some 
of the integrands near the integration boundary points. 
A midpoint scheme gives reasonable results, although, 
we found that an application of Aitken's method, see 
Ref. |31j . near the boundaries, combined with a trape- 
zoidal scheme in the non-divergent section yields more 
accurate and stable results. 



1. Ellipsoids 

In the case of an ellipsoidal particle there are three 
subdomains which partition D namely 



D n = [0,pi] x [0,tt/2]; 
D12 = \pi,P 2 ] x [0,7r/2]; 
D 2 = [P2,oo) x [0,tt/2], 



where 



Pi 

P2 



a cos i 



b 2 sin 2 * 



(B6) 
(B7) 
(B8) 

(B9) 
(BIO) 



These boundaries give 2det(0) an d a transition point in 
the integration domain, where there is a change in inte- 
gration kernel. It should be noted that the definition of 
subdomain here is slightly convoluted, since p\ and P2 
depend on <j>. The notation [0,pi] x [0,7r/2] means that 
for a specific <fi £ [0, tt/2) the z domain is the line segment 
[O,pi(0]._ 

Let us introduce the following parameters, which cor- 
respond to two coordinates of the plane-ellipsoid inter- 
section 



x± 



y± 



-b 2 z tan (j) ± ab\Jp\ — z 2 
(a 2 + b 2 tan 2 </>) cos <j> 



a 2 z ± a&tan <j>\Jp\ — ~Z' 
(a 2 + b 2 tan 2 0) cos 



(Bll) 



(B12) 



Using the above definitions, the semi-axes of the ellip- 
soidal cutout, and S12 are determined. The long semi- 
axis is given by 



(B13) 
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and the short semi-axis by 



b\ 1- 



li- 



2 a 



2b 



(B14) 



Let us further define the integral kernels 
Ii(rj) — 2nab. 



\ 



1 



1 




l\{rj) ( z-aq cos < 

J-2\n) — arccos 



(B15) 



(B16) 



b sin <j>\J\ — T) 2 
which can be applied to both oblate and prolate particles. 



Using Eqs. (B15 1 and (B16 1, the following expressions are 
obtained 



Ji 



J 2 



(v+/ a ) 
(v+/ a ) 

(»-/«) 



li(r))dr); 
X 2 (?7)d?7, 



(B17) 



(B18) 



which are related to the surface area Si, as we will show 
now. 

The total surface area of an ellipsoidal particle is now 
given by 



S' = / li(ri)dr], 



the surface area above the interface by 

Si(z,4>) = 



Ji + Ji (*, <f>) e D n 
J 2 (z, 4>) e d 12 
(z, 4>) e D 2 



and the cut-out surface area by 
Si 2 (z,(/>) ■■ 



7TQ cross ^cross 

o 



(z,<f>) G D 2 



(B19) 



(B20) 



(B21) 



To simplify the equation for the contact line length the 
following notations are introduced 



a 2 - b 2 

cross cross . 



(B22) 



-Across 

so that 
L(z,i 



= 4a r , 



tt/2 



l- C 2 ross sin 2 Vd^. (B23) 



Across (z, (/)) £ Dn U D12 

(z,<p)ED 2 



(B24) 



Equations (B22|-(B24| hold for both oblate and prolate 
ellipsoids, provided c 2 IOSS is allowed to assume negative 
values, when the particle is oblate. 



2. Cylinders 

For cylinders it is necessary to distinguish between two 
regimes in polar angle, separated by (f> — arctan(a/6), 
both of which have three z-domains. The angle <p gives 
the natural angle corresponding to the ratio of sides, 
which determines whether or not the plane can intersect 
the shaft of the cylinder without intersecting one of the 
end-caps. Some subdomains can be merged, which leads 
to the following partitioning of D 

Dn = [0,pi] x [0,0]; (B25) 

D13 - [0,pi] x [0,tt/2]; (B26) 

D 2 = \pi,p 2 ] x [0,7r/2]; (B27) 

D 3 = [^,00) x [0,tt/2], (B28) 

where 



Pi 



(B29) 



{a cos 4> — b sin cj) <f> € [0, 4>] 
b sin cj) — a cos <f> (j> £ [4>, it/ 2] 
p 2 = a cos <fi + 6 sin cf>, (B30) 



give the z-boundaries. These boundaries represent 
Zdct(0) and the position of the edge between the shaft 
and cap of a cylinder. 

The following parameters are introduced to aid nota- 
tion 



q+ 



z ± a cos i 
b sin <b 



(B31) 



u± = arccos q± — ?±\/l — q±, (B32) 



v± = q± arccos q± — y 1 — q±; (B33) 
w± = b 2 (u± ± 2u ± tan0). (B34) 

The parameter q± is related to the intersection of the 
plane and the end-cap, the others are parts of evaluated 
integrals. The integral kernel 



£«0 



2b 



1 -sin^sin 2 ^, (B35) 



is defined to help determine the contact line length. 

Using the above definitions, the following equations for 
S, Si, S12, and L, are derived. The total surface area is 



S = 2irb 2 + Airab. 



(B36) 



The area of the colloid's surface above the interface is 
given by 



Si{z,4>) = 



tt6 2 (1 - 2g_tan0) (z,4>) € Du 

w + + w- (z, (f>) e Di 2 

w- (z, 4>) e D 2 

(z,^)eD 3 



(B37) 



1G 



and the area of the cylinder-plane intersection is given 

by 



of equations" 



irb 2 

COS ( 



b 2 {u- — u + ) 



cos (, 
b 2 u_ 



(z,4>) € D u 



(z,4>) e D 12 



(B38) 



cos <p 
(z, 4>) G D 3 

The contact line length is found using the following "set 



L(z,< 



[ 1C(il>)d1> (z,4>)€D n 

J— it/2 

/>arcsinq+ ■ ■ 

/ IC(^)d^ + 2b^l-q 2 + + 2b^l-q 2 _ (z,(f>)eD 12 

J arcsin q- 

,7r/2 

/ /C(V)dV> + 2b^l - q 2 _ (z,4>)eD 2 

J arcsin o_ 



(B39) 



The above equations hold for both oblate and prolate 
cylindrical particles. 



3. Spherocylinders 

For spherocylinders the situation is even more compli- 
cated than it is for cylinders. Recall that in this paper 
we deviate from the classical definition of aspect ratio for 
spherocylinders, by including the caps in the length of a 
prolate particle and the toroidal rim in the width of an 
oblate particle. Hence, there are terms proportional to 
(a — b) present in the equations for prolate spherocylin- 
ders, which correspond to the traditionally used length, 
and (6 — a) terms for oblate particles, which correspond 
to the traditionally used width (also see main text). 

Since we did not succeed in formulating a single set 
of equations which holds for both oblate and prolate 
spherocylinders, we have split this paragraph into two 
parts. The first part describes prolate particles, the sec- 
ond oblate particles. For both species there are again two 
polar angle regimes, each having four z-rcgimcs. This 
makes the notation in the following quite heavy, espe- 
cially because a large number of parameters is introduced 
to formulate these equations as elegantly as possible. 



a. Prolate 

In the case of a prolate spherocylinder there are two po- 
lar angle regimes, separated by the angle <f> = arctan((a — 
b)/b), each of which can be split into four z-domains. 
Again </> is the natural transition angle related to the ratio 
of sides of the cylindrical part of the particle. After re- 
duction D can be written as D = _D 11 U-Di 2 U-D 2 U.D3UD4, 
with 





= [0,pi] x [0,0]; 


(B40) 


D 12 


= [0,pi] x [0,tt/2]; 


(B41) 


D 2 


= [pi.pa] x [0,7r/2]; 


(B42) 


D 3 


= [pa.ps] x [0,7r/2]; 


(B43) 


D± 


= [P3,oo) X [0,7r/2], 


(B44) 



where 



Pi 



{(a — b) cos — 6 sin <p <E [0, 0] 

fesin <f> — (a — b) cos0 G [0, it/ 2] 
p 2 = (a — b) cos0 + 6sin0; 
P3 = (a — b) cos (f> + b. 



(B45) 

(B46) 
(B47) 
(B48) 



The three boundary values for z correspond to z<iet(0)> 
the transition between the shaft and a partially inter- 
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sected sphere-cap, and the transition between to the sit- 
uation where only the sphere-cap is intersected and not 
the shaft. 

Let us now redefine some variables. Note that in the 
following q±, u±, and v± play a similar role as for the 
cylinder. 



s± 

t± 

u± 
v ± 
w± 
x± 

T 2 
X± 

y± 
a+ 

A_ 



z ± (a — b) cos < 
b sin d> 



cos cj)J 1 + (1 - q±) tan 2 <p ± q± sin 



arcsin q± + q± J 1 — gj_ 



arccos q± — q± \J 1 — q± ; 
(/± arccos g± 
<?± cos<Wl -<?|; 



1 — q 2 - sin 2 



1 - sin 2 . 



arccos 



g± cos < 
x± 



= X-(-K-y-); 

= x\y+; 

= x 2 _(iz-y-), 



(B49) 

(B50) 

(B51) 

(B52) 

(B53) 

(B54) 

(B55) 
(B56) 

(B57) 

(B58) 
(B59) 
(B60) 
(B61) 



where q± and s± are intersection related quantities, and 
the other quantities are to simplify the notation of evalu- 
ated integrals. The value x± is introduced here, because 
x± is defined as 1 — q\ sin 2 </>, rather than |1 — q\ sin 2 <j>\. 
The appearance of the absolute value in the latter would 
be a natural consequence of taking the square of x± . In 
this way the notation here is similar to c 2 ross for ellip- 
soids, for which we also used a "=" symbol. We realize 
that this notation is somewhat unconventional, but there 
are instances where x\ < in our definition. However, 
it should be pointed out that taking the square root to 
determine x± is not a problem, because on the domains 
that x\ is negative the equations do not contain instances 
of x± . The following integral kernels 



/c(V) 



2b 



1 — sin 2 </>sin 2 



£±{h) = 2b 2 arccos 



(h =F q±) tan < 



(B62) 
(B63) 



are (re)defined. Here the =F-sign in Eq. (B63 1 is the cause 
of an asymmetry which appears in the integration bound- 
aries later. 

Using the above equations and kernels the following 
physical quantities are derived, which determine the ad- 
sorption free energy 



S = Anab: 



(B64) 



Si(z,i 



27r6 2 (l — g_ tan ( 



= < 



C-{h)Ah- / C + (h)dh + nb%2-q + -q-) + 2b 2 (v + -V-)ta,nt 
[ £_(/i)d/i + 7r& 2 (l-q_)-26 2 t>_tan0 

J-q- 

2nb(p 3 - z) 




z,4>) e d 12 

z,4>) G D 2 

z,4>) e D 3 
z, (/)) e D A 



(B65) 
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b 2 (t+-t-) 
cos(0) 



7Tb 2 

cos(0) 

+ b 2 (^+ ~ w + ) + b 2 (fi- + w-) 



COS ( 



+ 6 2 (>_ + «;_) 



irb a; 



2^,2 



z,0) G D n 

z,4>) G D 12 

z, 0) G L> 2 
z, 0) G D 3 
z, 0) G Di 



(B66) 



L(z,( 



/•7r/2 



/ K(i>)dtp 

J-tt/2 


(z,</ 


) € D n 




/■arcsin g+ 

/ lC(ip)dip + 2b(\ + + \-) 

/arcsin q_ 


(z,</ 


) G Dia 




•/ arcsin q_ 


(z, 


) G £> 2 


(B67) 




(z,0 


) 6 ^3 







(z,r/ 


) G L> 4 





Note the sign asymmetry in the integration boundaries where 



of the C± integrals in Eq. ( B65 1 is induced by the =F-sign 
in Eq. ( |B63 1. Attempts to rewrite the integral in such 
a way that the asymmetry is eliminated lead to results 
which look contrived and are still asymmetric in a cer- 
tain way. Although the appearance of asymmetries may 
seem unphysical, we have extensively verified that these 
equations indeed hold. 



b. Oblate 

For oblate particles none of the integral equations, 
which describe the surface areas and contact line length, 
can be evaluated to obtain closed analytic expressions 
in terms of standard functions. Again 5 relevant do- 
mains are found, where the 0-domain is split by <p = 
arctan(a/(6 — a)). The angle is related to the dimen- 
sions of the cylindrical core of the prolate sphcrocylindcr. 
The subdomains are given by 



D u = 


[Q,pi] x [0,0]; 


(B68) 


D 12 = 


[Q,pi] x [0,tt/2]; 


(B69) 


D 2 = 


[ Pl ,p 2 ] x [0,tt/2]; 


(B70) 


D 3 = 


\P2,Pa] x [0,tt/2]; 


(B71) 


D± = 


[ P3 ,oo) x [0,tt/2], 


(B72) 



Pi = 

Vi = 

P3 = 



{a cos — (b — a) sin < 
(6 — a) sin — a cos < 
(6 — a) sin + a cos 0; 
(b — a) sin + a. 



0e [0,0] 

0G[0,7r/2] 



(B73) 

(B74) 
(B75) 
(B76) 



The three z-domain p.; indicate Zd e t(0)j and the position 
natural transition points on the particle's surface. 



To ease notation the following variables are introduced 
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4 
4 



cos 0(z — (b — a) sin 0) ± sin (f>\/ (a + z — (b — a) sin cj)) (a — z + (b — a) sin <p) ; 
cos </>(z + (b — a) sin 0) ± sin (a + z + (b — a) sin <fi) (a — z — {b — a) sin <p) ; 



^ = 27ra a - i ! ± + (6 - o) arccos 



(B77) 
(B78) 

(B79) 



where i can be cither T or '6'. These stand for 'top' and 
'bottom' respectively, but the latter is in no way related 
to the rotational radius of the particle, which is also given 
by b. The points x*± and x b ± are locations where the 
interfaces intersects the spherocylinder in a conveniently 
chosen coordinate frame. The quantities [i % ± originate 
form the evaluation of integrals. The following useful 
functions arc defined to aid notation 



x — z cos ( 



z sin < 



tan < 



h(x) 

k(x) = (b — a) 2 arccos 



h(x) 



r(x) 
w(x) 



b — a , 

-h(x) a) 2 -h 2 (x). 

= (b — a) + \J a 2 — x 2 
= yjr 2 {x)-h 2 (x) 



(B80) 



(B81) 

(B82) 
(B83) 



Here h stands for a height related function, k is a inte- 



gration reduction function, r is a radial distance func- 
tion, and w is a width function. In this paragraph 
a prime denotes a derivative with respect to x, e.g., 
h'(x) = dh(x)/dx. For the oblate spherocylinder three 
integral kernels are required, namely 



2ar(x) 



Va 2 ^ 
2w(x) 



arccos 



h(x) 



r(x) 



M 2 {x) 

Mz{x) = 2y/l + (h'(x)) 2 + (w'(x)) 



sm< 



(B84) 

(B85) 
(B86) 



The first of the integral kernels is used in determining 
the area Si, the second to determine Si 2 , and the third 
is used to determine the contact line length. 

The above equations and kernels arc applied to derive 
the following equations 



S = 2tt (b 2 + (tt - 2)ba - (tt - 3)a 2 ) ; 



(B87) 



Si(z,, 



I N\(x)&x + Ti(b - a) 2 + 4 

J X 1 

WiO)dx + k(a) + fc(-a) 

a 

Afi(x)dx + fe(a) 

/ A^(x)dx 




z,<j>) g Dn 
0, 0) g z? 2 

0) G As 
0, <j)) G £>4 



(B88) 
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S 12 (z, 



= < 



Af 2 (x)dx (z,(j>) e D n 

f N 2 {x)dx {z,<j>)eD 12 

N 2 {x)dx (z,</>)eL> 2 

Af 2 (x)dx {z,4>) e D 3 

(z, 4>) e L>4 



(B89) 



L(z,</>) = < 



Af 3 (x)dx (z,<j>)eD n 

Af 3 (x)dx + 2w(a) +2w(-a) (z,<j>) £ D 12 

f Af 3 (x)dx + 2w(a) {z,cj>)£D 2 ' 

r x + 

/ AA 3 (x)dx (z»e£ 3 

(z,c/))eD 4 



(B90) 



r 



Note that although the above equations are more sym- still a degree of asymmetry in the boundary conditions, 
metric than those for prolate spherocylinders, there is 
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